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We study the formation of large vortex aggregates in a rapidly rotating dilute-gas Bose-Einstein 
condensate. When we remove atoms from the rotating condensate with a tightly focused, resonant 
laser, the density can be locally suppressed, while fast circulation of a ring-shaped superflow around 
the area of suppressed density is maintained. Thus a giant vortex core comprising 7 to 60 phase 
singularities is formed. The giant core is only metastable, and it will refill with distinguishable single 
vortices after many rotation cycles. The surprisingly long lifetime of the core can be attributed to 
the influence of strong Coriolis forces in the condensate. In addition we have been able to follow the 
precession of off-center giant vortices for more than 20 cycles. 
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The phenomenon of quantized vortex formation is a 
unifying feature present in many quantum mechanical 
systems. Vorticity is intimately connected to superflu- 
idity and thus is a premier means for fundamental and 
comparative studies of different "super" systems. In di- 
lute gas Bose-Einstein condensates, vortex experiments 
have ranged from the study of individual or few vortices 
[1-4] and vortex rings [5,6] to the first creation of vor- 
tex lattices [7] and the study of systems containing large 
amounts of vorticity [8,9]. 

Using a refinement of our experimental technique, we 
have now been able to generate a giant vortex. We will 
use the term "giant vortex" to denote a region, con- 
taining multiple units of vorticity, in which the density 
is completely suppressed such that the individual vor- 
tices are no longer discernible. We note that some au- 
thors reserve this term for higher order phase singulari- 
ties. The possible existence of stable giant vortices un- 
der certain conditions like supersonic flow in potentials 
that have strong quartic terms has been theoretically pre- 
dicted (Ref. [10-12], see also Ref. [13]), and in Ref. [14] 
it is shown that a pinning potential can lead to stable 
multiquantum vortices. In our case, however, the giant 
vortex formation arises as a dynamic effect. Nevertheless 
the lifetime of our giant vortices can extend over many 
seconds, which we attribute to a stabilization of density 
features in a rapidly rotating condensate due to strong 
Coriolis forces. The influence of Coriolis forces can also 
induce oscillations of the giant vortex core size in the 
early stages of its evolution. 

Our experimental starting conditions are similar to 
the ones described in our previous papers [9,15]. We 
start with a cloud of uncondensed 87 Rb atoms in the 
\F = l,mp = —1) state, held in a magnetic TOP trap 
and precooled to a temperature slightly above the critical 
temperature for Bose-Einstein condensation. By ellipti- 
cally deforming the magnetic trap in the horizontal xy 



plane and suddenly changing the angle of the deforma- 
tion, we impart angular momentum around the vertical 
z-axis to the cloud. We then restore axial symmetry and 
perform a quasi one-dimensional evaporation in a trap 
with frequencies {td p ,uj z } = 27r{8.3, 5.4}Hz. This evap- 
oration preferentially removes atoms close to the axis of 
rotation (z-axis), thus spinning up the cloud. At the end 
of this step, we routinely achieve fast rotating conden- 
sates consisting of 3 x 10 6 atoms with a rotation rate of 
95% of the radial trap frequency, with no detectable ther- 
mal cloud. This rotation rate is deduced from the aspect 
ratio, which changes from the value of 1.53 for a static 
cloud to 0.48 because of the centrifugal forces created by 
the rotation. These condensates typically contain 180 or 
more vortices as seen in Fig. 1(a) or Fig. 2(a). After care- 
ful optimization of the trap roundness and in the pres- 
ence of a quasi 1-D rf-shield we are able to observe the 
BEC rotation for times exceeding 5 minutes. While we 
do lose atoms from the condensate over this time scale, 
the rotation rate remains at its initial value. 

The evaporative spin-up mechanism is limited because 
eventually all the uncondensed atoms have been removed 
or condensed during the quasi one-dimensional evapora- 
tion. However, the angular momentum per particle in 
the condensate can be further enhanced by selectively 
removing atoms with low angular momenta by a reso- 
nant, focused laser beam sent through the condensate on 
the axis of rotation. Since the rotating BEC typically has 
a Thomas-Fermi radius of 66/im, the width of the laser 
beam of about 16/xm, stated here as the full width at half 
maximum (FWHM) of the Gaussian intensity profile, is 
small enough to provide sufficient selectivity. The fre- 
quency of the laser is tuned to the F" = 1 — * F' = tran- 
sition of the D2 line, and the recoil from spontaneously 
scattered photons blasts atoms out of the condensate. 

Using this new technique we are able to substantially 
increase the vorticity in the BEC and can now routinely 
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obtain condensates containing more than 250 vortices 
[Fig. 1(b)]. Since these vortex cores are too small to 
be imaged when the BEC is held in trap, the images in 
Figs. 1, 2, 3, and 6 are taken after having let the BEC 
expand. For this, a rapid adiabatic passage technique 
is employed to change the hyperfine state of the atoms 
from the original trapped \F = 1, nip = — 1) state to the 
antitrapped \F = 2 1 mp = —1) state. Subsequently the 
minimum of the magnetic trapping field is rapidly moved 
below the position of the atoms so that the atoms are sup- 
ported against gravity while at the same time they are 
radially expelled [16]. After a chosen expansion time the 
magnetic fields are switched off and the atoms are imaged 
along the original axis of rotation. This technique allows 
us to expand the condensate to an adjustable diameter 
that can exceed 1.4 mm [17]. 

If the removal of central atoms is driven strongly 
enough, a giant vortex appears. The core is surrounded 
by a ringshaped supcrflow. The circulation of this su- 
pcrflow is given by summing up all the "missing" phase 
singularities in the core and can assume very large values. 
Using the fact that for large amounts of circulation the 
rotation of the velocity field can be obtained classically, 
we determine that in some cases the BEC supports a ring- 
shaped superflow with a circulation of up to 60 quanta 
around the core. We can easily control the amount of 
this circulation by changing the intensity or the duration 
of the laser beam that removes atoms along the axis of 
rotation. 

In our experiment the formation of the giant vortex 
comes about in a sequence of very distinct stages as 
shown in Fig. 1. For this expansion image sequence a 
rapidly rotating BEC is first formed by our evaporative 
spin-up technique. Then the atom-removal laser is ap- 
plied with a fixed power of 8 fW for a variable amount 
of time, followed by a 10 ms in-trap evolution time and a 
45 ms expansion in the antitrapping configuration de- 
scribed above. During the atom removal an rf-shicld 
is left on to constantly remove uncondensed atoms that 
tend to decelerate the condensate rotation. Figure 1(a) 
shows the result of only the evaporative spin-up. This 
particular condensate contains 180 vortices and has a 
Thomas-Fermi radius of 63.5 fxm when held in the trap. 
When the atom removal laser is applied for 14 s as in 
Fig. 1(b), the number of vortices is increased to 250 
and the Thomas- Fermi radius to 71 \xm. The rotation 
rate, determined from sideview aspect-ratio images, has 
increased from 0.94 uj p to 0.97 uj p . After atom removal 
times of 15 to 20 s, the vortex lattice becomes disor- 
dered [Fig. l(c,d)], and the giant vortex core starts to 
develop in the center [Fig. l(e,f)]. An enlarged view of 
the core region in this early stage of giant vortex core 
formation is seen in Fig. l(i), which nicely demonstrates 
how the individual vortices in the center consolidate. For 
the larger removal times shown in Figs. 1(g) and (h), a 
clear elliptical deformation is observed. This deforma- 



tion is due to residual static asymmetries of the trapping 
potential that can nearly resonantly drive quadrupolar 
surface modes when the rotation frequency of the con- 
densate approaches the radial trapping frequency. We 
have verified this explanation by deliberately changing 
the roundness of the trapping potential and observing 
corresponding changes in the ellipticity and in the direc- 
tion of the ellipse. If the trap roundness is not optimized, 
the ring-shaped supcrflow has a tendency to fragment 
into typically 5 to 8 blobs that continue to rotate around 
the core. 

If the atom removal laser is applied only for a limited 
amount of time, the giant vortex will fill in and the vor- 
tices rcfrccze [18] into a well ordered lattice again after 
a sufficiently long evolution time. An example is shown 
in Fig. 2. Here, the starting point is the condensate in 
Fig. 2(a), containing 190 vortices and rotating at 0.95 uj p . 
A pulse from the atom removal laser creates the core 
seen in Fig. 2(b). Crudely assuming that all atoms, and 
only those atoms, originally within a cylindrical volume 
of one-third of the original cloud radius are removed, we 
can perform an integral over the Thomas-Fermi density 
profile with a rigid-body-rotation velocity distribution, 
and find that 23% of the atoms have been removed but 
only 4% of the angular momentum. Since in equilibrium 
the total number of vortices should be linear in the aver- 
age angular momentum per atom [19], we would expect 
that the original 190- vortex cloud should reequilibrate to 
form a cloud with 236 vortices, in rough agreement with 
the value of 260 vortices observed in the reequilibrated 
cloud shown in Fig. 2(d). 

In our case the formation of the giant vortex is not 
connected to a repulsive conservative optical potential 
produced by the central laser beam, but only due to the 
removal of atoms from the axis of rotation by sponta- 
neous photon scattering. To prove this, we have var- 
ied the laser frequency over the range from -6 MHz to 
+6 MHz around the F" = 1 — > F' = resonance and in 
all cases have been able to generate a giant vortex. 

Since theoretical studies predict the appearance of a 
giant vortex if the trapping potential has a sufficiently 
strong admixture of quartic terms, it is important to 
quantify the strength of the quartic terms in our case. 
For this we have calculated the term of the TOP trap 
potential that depends on the radial position as r 4 . At 
the outer edge of the condensate this quartic term pro- 
duces a correction to the potential of less than 10 -3 , com- 
pared to the effective r 2 term, even though the latter is 
much weakened by the centrifugal force. Higher order 
corrections contribute even less. Empirically, the refill- 
ing and refreezing we observe (Fig. 2) indicates that the 
giant vortex is not an equilibrium configuration but has 
instead only dynamical stability. 

We suggest an intuitive, classical picture describing the 
formation as a dynamical effect: the removal of atoms 
from the center of the condensate produces a pressure 
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gradient due to mean field energy that tries to drive 
atoms from the outer regions into the center, so as to 
close the hole. Due to Coriolis forces, however, atoms 
moving radially towards the center are deflected and as- 
sume a fast azimuthal motion around the core rather than 
filling the core, thus creating the giant vortex. 

Some exotic features of the early stages of core forma- 
tion are revealed when we apply only a very short, weak 
atom removal pulse and observe the subsequent evolu- 
tion as in Fig. 3. Here, the atom removal laser has a 
power of 2.5 pW and a fixed pulse length of 5 ms. Its 
FWHM of 16 /um is approximately twice the lattice spac- 
ing, 7/zm. By varying the delay time between the laser 
pulse and the expansion, we see (Fig. 3) that the core 
formation clearly lags behind the laser pulse. Fig. 3(e) 
shows a zoomed- in view of the core region of Figure 3(c). 
It is very interesting to observe that for these short, weak 
pulses the density depression appears to develop in dis- 
crete steps, and the step boundaries follow the hexagonal 
lattice pattern. 

When applying stronger laser pulses, we detect clear, 
damped oscillations of the core area as shown in Fig. 4. 
We study these oscillations by analyzing in-trap images 
taken after a variable evolution time between the end 
of the atom-removal pulse and the expansion, such as 
the ones shown in the inset of Fig. 4. The oscillation 
frequency depends on the initial conditions; for the case 
of Fig. 4 with initially 2.2 x 10 6 atoms at a rotation rate 
of 0.9 w p we obtain a frequency of 3.5 u> p . Decreasing the 
initial rotation rate of the condensate leads to faster core- 
oscillation frequencies and increases the amplitude of this 
oscillation; for rotation rates below 0.9 lo p in-trap images 
even show a near complete closure of the core, followed by 
the core opening again (see inset of Fig. 4). Presumably 
these oscillations arise when a sudden removal of atoms 
leaves forces from the density gradient and Coriolis forces 
initially out of equilibrium. The core oscillation is not 
observed to be related to an overall breathing mode of 
the condensate. 

In addition we can observe the precession [2] of an off- 
center giant vortex. For this, the atom removal laser 
is deliberately offset from the center of cloud rotation. 
The duration of the removal pulse is kept short (10 ms) 
in comparison to the initial vortex lattice rotation pe- 
riod (27r/(0.95 Wp) = 126 ms) so as to create only a local 
hole. By using different laser powers, we can vary the 
size of the core from a small hole as in Fig. 5(a) all the 
way to the extreme case of a big hole that only leaves 
a crescent segment of the BEC as shown in Fig. 5(c), 
and the precession of these holes can be monitored by 
applying variable evolution times of the trapped BEC af- 
ter the atom removal pulse is finished. In all three cases 
the measured precession frequency is approximately ui p . 
In the case of Fig. 5(b) and (c) we are able to follow 
the precession for more than 20 cycles. If, instead, the 
laser is left on for approximately a full lattice rotation 



cycle, a complete ring can be cut out of the condensate 
[Fig. 6(a)], which eventually breaks up into many indi- 
vidual blobs [Fig. 6(c)], each of which supports remnants 
of the original lattice. Both the long-lived core preces- 
sion as well as the ring structure in Fig. 6 are impressive 
demonstrations of the stability of density features due to 
Coriolis forces in rapidly rotating condensates. 

Finally, we also observe that the sudden extraction of 
atoms can launch transverse modes in the lattice [20], 
which will be the subject of a future publication. 

The work presented in this paper was funded by NSF 
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FIG. 1. Different stages of giant vortex formation process, 
(a) Starting point: BEC after evaporative spin-up. (b)-(h) 
Laser shone onto BEC for (b) 14 s, (c) 15 s, (d) 20 s, (e) 22 s, 
(f) 23 s, (g) 40 s, (h) 70 s. Pictures are taken after 5.7-fold 
expansion of the BEC. Laser power is 8 fW. (i) Zoomed-in 
core region of (f). 
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FIG. 3. Core developing after a 5 ms short, 2.5 pW laser 
pulse. In-trap evolution time after end of pulse is (a) 0.5 ms, 
(b) 10.5 ms, (c) 20.5 ms, (d) 30.5 ms. Images taken after 
sixfold expansion of BEC. (e) Zoomed-in core region of (c). 
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FIG. 2. Lattice reforming after giant vortex formation, (a) 
BEC after evaporative spin-up. (b) Effect of a 60 fW, 2.5 s 
laser pulse. (c),(d) Same as (b), but additional 10 s (c) and 
20 s (d) in-trap evolution time after end of laser pulse. Images 
taken after a sixfold expansion of the BEC. 



FIG. 4. Oscillation of core area after an 8 pW, 5 ms laser 
pulse. Starting conditions are 2.5 x 10 atoms with rotation 
rate 0.9 uj p . Time given is in-trap evolution time after end of 
pulse. Core area is normalized to mean of all data points. 
Inset: Initial conditions 3.5 x 10 6 atoms with rotation rate 
0.78 ui p ; in-trap images taken after a 14 pW, 5 ms laser pulse 
followed by evolution time of (a) 20 ms, (b) 40 ms, (c) 60 ms, 
and (d) 80 ms. 



4 



55 ms 



85 ms 



115 ms 145 ms 1125 ms 3625 ms 



a) 


9 


4» 


i 


□ 




* 


b) 


T 


B 


c 


D 




□ 


% 


> 






H 







FIG. 5. Giant core precession. Cores created by a 10 ms 
off-centered laser pulse with a power of (a) 4.2 pW, (b) 33 pW, 
(c) 470 pW. Time given in figure is in-trap evolution time 
after end of laser pulse. For reference, a naive expectation for 
the lifetime of a giant vortex given by the radius of the core 
divided by the speed of sound in the surrounding cloud would 
be 13 ms and 30 ms for (a) and (b), respectively. 




FIG. 6. Ring cut out of a BEC by a 125 ms long, 
off-centered laser pulse. Expansion image taken (a) directly 
after end of pulse (b) after an additional in-trap evolution 
time of 200 ms after end of pulse and (c) evolution time 2 s. 
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